arXiv:1505.02881vl [astro-ph.GA] 12 May 2015 


Publ. Astron. Soc. Japan {20^4) 00(0), 1-5 
doi: 10.1093/pasj/xxx000 


1 


Radial distributions of arm-gas offsets as an 
observational test of spiral theories 

Junichi Baba\ Kana Morokuma-Matsui^, and Fumi Egusa^^ 

^ Earth-Life Science institute, Tokyo institute of Technoiogy, Ookayama, Meguro, Tokyo 
152-8551, Japan. 

^ Chiie Observatory, Nationai Astronomicai Observatory of Japan, 2-21-1 Osawa, Mitaka, 
Tokyo 181-8588, Japan. 

^ institute of Space and Astronauticai Science, Japan Aerospace Expioration Agency, 
Sagamihara, Kanagawa 252-5210, Japan 

*E-maii: babajn@eisi.jp; babajn2000@gmaii.com 
Received <09-Apr-2015>; Accepted {12-May-2015> 


Abstract 

Theories of steiiar spirai arms in disk gaiaxies can be grouped into two ciasses based on the 
iongevity of a spirai arm. Aithough the quasi-stationary density wave theory supposes that 
spirais are rigidiy-rotating, iong-iived patterns, the dynamic spirai theory predicts that spirais 
are differentiaiiy-rotating, transient, recurrent patterns, in order to distinguish between the 
two spirai modeis from observations, we performed hydrodynamic simuiations with steady and 
dynamic spirai modeis. Hydrodynamics simuiations in steady spirai modeis demonstrated that 
the dust iane iocations reiative to the steiiar spirai arms (hereafter, arm-gas offsets) depend on 
radius, regardiess of the strength and pitch angie of the spirai and the modei of the inter-steiiar 
medium (iSM). in contrast, we found that the dynamic spirai modeis show no systematic radiai 
dependence of the arm-gas offsets. The arm-gas offset radiai profiie method, together with the 
other test methods, wiii heip us to distinguish between the two spirai modeis in observed spirai 
gaiaxies. 

Key words: Gaiaxies; spirai — Gaiaxies: kinematics and dynamics — ISM: kinematics and dynamics — 
method: numerical 


1 Introduction 

Spiral structures in galaxies have been interpreted as “standing 
waves” in a stellar disk (the so-called quasi-stationary density 
wave hypothesis; Lin & Shu 1964; Beitin & Lin 1996). The 
density waves generate spiral perturbations and affect on the 
gas flows. When the gas inside a co-rotation radius overtakes 
a spiral arm as it moves around a galactic disk, the gas is ex¬ 
pected to form a standing shock, called a “galactic shock”, at the 
upstream side of the spiral arm (e.g., Fujimoto 1968; Roberts 
1969; Shu et al. 1973). Since such a galactic shock strongly 
compresses the gas, the galactic spiral dust-lanes and the as¬ 
sociated star formation around spiral arms have since been re¬ 
garded as consequences of the galactic shocks. The galactic 


shock hypothesis qualitatively explains the CO-Ha offset dis¬ 
tributions observed in some spiral galaxies such as M51 (e.g., 
Egusa et al. 2009; Louie et al. 2013). However, it is still un¬ 
clear whether a spiral arm rotates with a single pattern speed 
for several rotational periods (i.e., ~ 1 Gyr), as predicted by the 
quasi-stationary density wave hypothesis. 

On the other hand, a transient spiral hypothesis was pro¬ 
posed in the 1960s (Goldreich & Lynden-Bell 1965; Julian & 
Toomre 1966). This hypothesis has been extended to a tran¬ 
sient recurrent spiral (hereafter, “dynamic spiral”) hypothesis 
by iV-body simulations of stellar disks (see a review by Dobbs 
& Baba 2014); the amplitudes of stellar spiral arms change 
on the time scale of galactic rotation or even less (i.e, a few 
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hundreds of Myrs; Sellwood & Carlberg 1984; Baba et al. 
2009; Fuji! et al. 2011; Sellwood 2011; Grand et al. 2012; Baba 
et al. 2013; D’Onghia et al. 2013; Pettitt et al. 2015), and 
the spiral at any given radius co-rotates with materials (Wada 
et al. 2011; Grand et al. 2012; Baba et al. 2013; Kawata et al. 
2014). Fuji! et al. (2011) revealed a self-regulating mechanism 
that maintains multi-arm spiral features for at least 10 Gyrs. 
The “dynamic equilibrium” nature of stellar spiral arms can 
be attributed to their co-rotational nature (Baba et al. 2013). 
Furthermore, Dobbs & Bonnell (2008) and Wada, Baba & 
Saitoh (2011) performed hydrodynamic simulations of the gas 
flows in dynamic spirals, and showed that gas does not flow 
through a spiral arm, but rather effectively falls into the spiral 
potential minimum from both sides of the arm (called “large- 
scale colliding flows”). The dynamic spiral hypothesis is also 
supported by observations (e.g., Foyle et al. 2011; Ferreras et al. 
2012 ). 

Some observational tests for discriminating between quasi¬ 
stationary density waves (i.e., steady spirals) and dynamic spi¬ 
rals have been proposed (see a review by Dobbs & Baba 2014). 
In this paper, we propose a new indicator to discriminate be¬ 
tween the two spiral models based on the radial distributions of 
arm-gas offset angles. The dust lanes locations relative to a stel¬ 
lar spiral arm depend on the gas flow in the spiral galaxy. In the 
following sections, we demonstrate that the arm-gas offset dis¬ 
tributions show different radial dependences between hydrody¬ 
namic simulations of steady spiral and dynamic spiral models. 
This difference suggests that measuring the radial distribution 
of the arm-gas offset angles can be a new observational means 
of distinguishing between the two spiral models. 


2 Models and methods 

We performed hydrodynamic simulations in rigidly rotating spi¬ 
ral potentials and W-body/hydrodynamic simulations of stel¬ 
lar and gas disks with the W-body/smoothed particle hydro¬ 
dynamics (SPH) simulation code ASURA-2 (Saitoh & Makino 
2009,2010). 

2.1 Steady spiral models 

To simulate the gas flows in spiral galaxies with stationary 
density waves, we impose the rigidly rotating spiral potentials 
into the axisymmetric potential. The gravitational potential 
of the model galaxy consists of an axisymmetric and a non- 
axisymmetric (i.e., spiral) parts. The static axisymmetric part of 
the gravitational potential comprises of a stellar disk, a spher¬ 
ical bulge and halo. The initial radial distribution of the gas 
follows an exponential profile with a scale length of ~ 11 kpc, 
which is motivated by the observations of the Milky Way galaxy 
(e.g., Bigiel & Blitz 2012). Figure 1 shows the circular velocity 


curves of each component. 

The gravitational potential of the spiral arm is given by 


= A{R,z) 
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where A, m, isp, Dp, and zo are the amplitude of the spiral po¬ 
tential, the number of stellar spiral arms, the pitch angle, the 
pattern speed, and the scale-height, respectively. Note that 4?sp 
decreases from R=1 kpc exponentially to be <l?sp = 0 at the cen¬ 
ter. In this study, we set m = 2, zq = 100 pc, and i?o = 1 kpc. 
We adopt a typical value of the pattern speeds of spiral arms 
in Milky Way galaxy. Dp ~ 23 km s“^ kpc“^ (i.e., Ren = 10 
kpc), although they have a wide range of 18 - 30 km s“^ kpc“^ 
(Gerhard 2011). The amplitude of the spiral potential is con¬ 
trolled by the dimensionless parameter 


, ( 2 ) 

|®o|sintsp 

where T'o(i?,2) is the axisymmetric potential. This dimension¬ 
less parameter measures the gravitational force due to the spiral 
arms in the direction perpendicular to the arms relative to the 
radial force from the background axisymmetric potential (Shu 
etal. 1973). 

In the steady spiral runs, we perform two series of runs 
with the different ISM models: the first series is hydrodynamic 
simulations of an isothermal, non-self-gravitating gas. No star 
formation and stellar feedback are considered. The temper¬ 
ature of the gas is assumed to 10'* K to effectively include 
a contribution of turbulent motions. This ISM model corre¬ 
sponds to that assumed in the classical galactic shock theory 
(e.g., Roberts 1969). The second series is hydrodynamic sim¬ 
ulations with a more realistic ISM model, which includes the 
self-gravity, radiative cooling for a wide temperature range of 
20 K < T < 10® K (Wada et al. 2009), heating due to far- 
ultraviolet radiation (Gerritsen & Icke 1997), star formation and 
supernova feedback (Saitoh et al. 2008). The self-gravity is 
calculated by the Tree with GRAPE method with a software 
emulator of GRAPE, Phantom-GRAPE (Tanikawa et al. 2013). 
Model parameters are summarized in Table 1. The simulations 
are performed in a frame co-rotating with the spiral arms until 
t = 350 Myr. To avoid strong transients in the gas flows caused 
by a sudden introduction of the spiral potential, we increase its 
amplitude linearly from t — 100 Myr to 200 Myr. 


2.2 Dynamic spiral models 

We performed three-dimensional, high-resolution M-body/SPH 
simulations of a Milky Way-like galaxy model. An initial ax¬ 
isymmetric model of stellar and gaseous disks, a classical bulge 
embedded in a dark matter halo has the circular velocity curves 
of each component shown in Figure 1. We take into account 
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Fig. 1. Circular velocity curves of each component. In the cases of the 
rigidly rotating spiral arms, all components are assumed to be external 
potentials (section 3.1). For the case of the dynamically evolving 
barred-spiral, only the dark matter halo is treated as an external potential, 
but other components are solved by the Af-body/SPH methods (section 3.2) 

the self-gravity of the gas, cooling, heating, star formation 
and supernova feedback (section 2.1). In this model, the stel¬ 
lar bar is developed by a bar instability (Ostriker & Peebles 
1973; Efstathiou et al. 1982) at f ~ 1.5 Gyr. In this paper, we 
focus on the gas flows and spatial distributions in simulated un¬ 
barred (t ~ 1 Gyr) and barred {t ~ 2.5 Gyr) spiral galaxeis. 
Detailed studies on dynamics of stars, as well as structures of 
the ISM/clouds, will be separately discussed in forthcoming pa¬ 
pers (J. Baba et al., in preparation). 

Table 1 . Model parameters for the steady spiral models and prop¬ 
erties of spiral arms in dynamic spiral models. Note that the ‘mul¬ 
tiphase’ ISM model includes cooling, heating, star formation and 
supernova feedback (See section 2 for details). For comparison, 
typical values of iap and J' for dynamic spiral models are pre¬ 
sented. 


Models 

^sp 

[deg.] 

T 

[%] 

gas self 
-gravity 

ISM 

model 

steady spirals 

Sil0F02iso 

10 

2 

no 

isothermal 

SilOFOSiso 

10 

5 

no 

isothermal 

Si20F05iso 

20 

5 

no 

isothermal 

Si20F05msg 

20 

5 

yes 

multiphase 

dynamic spirals 

DynUnBar 

~25 

<5 

yes 

multiphase 

DynBar 

-30 

<8 

yes 

multiphase 


3 Results 

3.1 Arm-gas offset angles in steady spirals 

Our results show that the locations of galactic shocks strongly 
depend on radius. Figure 2 shows the radial distributions of 
the offset angle of the dense gas relative to the spiral potential 
(hereafter, the arm-gas offsets). In the tightly winding, weak 


spiral case (Sil0F02iso model; panel A of Figure 2), the galactic 
shock is located on the upstream side except for R<2 kpc. This 
is consistent with the predictions by the classical galactic shock 
theory (Roberts 1969). However, the galactic shock is away 
from the potential minimum and moves upstream if it is near 
the co-rotation radius (i?cr = 10 kpc). This is because the Mach 
number is a function of radius (e.g., Shu et al. 1973). Such 
radial dependence is consistent with Gittins & Clarke (2004), 
who found that the pitch angle of gaseous arms is smaller than 
that of stellar arms (see also Kim & Kim 2014). 

The radial trend of arm-gas offset angles in steady spirals 
depends qualitatively on neither the strength nor the pitch angle 
of the spiral. If we compare results with the different T mod¬ 
els (panels A and B in Figure 2), the galactic shock locations 
in SilOFOSiso model are shifted farther downstream (see also 
Woodward 1975). Solid green lines in panels B and C in Figure 
2 indicate the locations of galactic shocks in the Sil0F05iso 
model. The locations of galactic shocks in both models are al¬ 
most the same. This is a natural consequence of the definition 
of T. Note that the shocked layers in spiral potentials are not 
always dynamically stable (so-called wiggle instability). This is 
consistent with the previous time-dependent, multi-dimensional 
hydrodynamic simulations (e.g., Wada & Koda 2004; Shetty & 
Ostriker 2006; Dobbs & Bonnell 2006; Kim et al. 2014). The 
physical origin of wiggle instabilities is a standing problem', 
although the wiggle instabilities result in forming the substruc¬ 
tures such as spurs in the inter-arm regions. 

The steady spiral models with more realistic ISM models 
also show the similar radially dependent arm-gas offset angles. 
In panel (D) of Figure 2, we show the radial distribution of 
the arm-gas offsets for Si20F05msg model, which include self¬ 
gravity of the gas, radiative cooling, heating and stellar feed¬ 
back. As suggested by Lee (2014), the self-gravity of the gas 
may shift the shock location downstream (Sawa 1977; Wada 
2008), but a clear radial dependence of the arm-gas offset is 
still seen in Si20F05msg model. This suggests that an arm-gas 
offset distribution does not depend on ISM modeling, but is es¬ 
sentially determined by dynamical properties of a spiral arm. 

3.2 Arm-gas offset angles In dynamic spirals 

We investigate the dust-lane locations in dynamically evolving 
spiral arms. Figure 3 shows the B-band brightness distributions 
of stars with dust extinction at f = 1.12 Gyr and t = 2.55 Gyr. 
The model galaxy is an ‘unbarred’ grand-design spiral at f < 1.5 

' Hanawa & Kikuchi (2012) pointed out that the wiggle instability (Wl) may 
be due entirely to numerical artifacts arising from the inability to resolve a 
shock inclined to numerical grids. However, Kim et al. (2014) found that 
the Wl is a physical phenomenon and originating from the generation of 
potential vorticity at a deformed shock front, rather than Kelvin-Helmholtz 
instabilities as proposed by Wada & Koda (2004). We do not aim to inves¬ 
tigate the physical origins here. See these papers for more details. 
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Fig. 2. Panels (A)-(D): The azimuth angle of the gas relative to the potential minimum (A<^sp) as functions of radius for models with steady spirals at t = 350 
Myr. The potential minimum of each radius is located at A0sp = o. The negative value of an offset angle (A0sp < o) corresponds to a leading side (i.e., 
downstream inside the co-rotating radius). The dotted (black) lines indicate the locations of shocks in the Si10F02iso model (panel A). The solid (green) lines 
indicate the locations of shocks in the SilOFOSiso model (panel B). Panels (E) and (F): Same as the panels (A)-(D), but for dynamic spiral models at f = 1.12 
Gyr (panel E) and 2.55 Gyr (panel F). The model galaxy is an 'unbarred' grand-design spiral at t < 1.5 Gyr, but is a 'barred' grand-design spiral after then (J. 
Baba et al. in preparation). The locations of the grand-design spirals are defined by the phases of the dominant Fourier modes (i.e., m = 2 or 3). The 
background gray-scale map presents a stellar density. The grand-design spiral is located at A(>Bp = o. 


Gyr, and then becomes a ‘barred’ grand-design spiral. At t = 
2.55 Gyr, a stellar bar is seen at R<3 kpc. Properties of spirals 
are analyzed by using the ID Fourier decomposition method as 
follows: 

Am{R)^ J (3) 

where p{R, (j)) is the stellar mass distribution and p{R) is the 
azimuthally averaged surface stellar density at radius R, and 
Am{R) denotes a complex amplitude of the m-th mode, respec¬ 
tively. We analyze angular positions and angular phase speeds 
of spiral arms with the Fourier density peak method (Wada et al. 
2011; Roca-Fabrega et al. 2013). 

The angular phase speeds of spiral arms in dynamic spiral 
models follow the galactic rotation at almost all radii or radi¬ 
ally decrease (bottom panels of Figure 3). This is similar to the 
simulated multi-armed spirals (e.g.. Baba et al. 2013), as well 
as the simulated barred spirals (e.g.. Grand et al. 2012), which 
have the radially dependent rotation of spiral arms, but slightly 
faster than the galactic rotation. This suggests that spiral arms 
even in barred galaxies could not be rigid-body rotating patterns 
predicted by the quasi-stationary density wave theory. 

In contrast to the steady spiral models, we can clearly see 
that the dust lanes are along with stellar spiral arms in the dy¬ 
namic spiral model (top panels of Figure 3). As shown in the 
panels (E) and (F) of Figure 2, the radial distributions of the 
arm-gas offset angle show no clear radial dependence in dy¬ 
namic spiral models. This is because the gas does not flow 
through the spiral arm, and the gaseous arm remains until the 
stellar arm disperses (see also Dobbs & Bonnell 2008; Wada 
et al. 2011). 
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Fig. 3. Upper: B-band brightness distributions of stars (gray) in dynamic 
spiral models at t = 1.12 Gyr (left) and 2.55 Gyr (right). The dust extinction 
for the B-band map is estimated by multiplying a factor of , where the 
optical depth tb is calculated from the absorption cross section 
CTB = 6 X 10“^^ cm^ and the total hydrogen column density Eh. Lower: 
Radial profiles of the Fourier amplitudes and angular phase speeds of the 
dominant modes (m = 3 and 2) at t = 1.12 Gyr (left) and 2.55 Gyr (right), 
respectively. Solid line indicates a galactic rotation frequency n, and 
dashed lines indicate n ± n/m, where k . is an epicyclic frequency. The 
angular speeds are not well determined around radii with weak amplitudes. 
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4 Discussion and Summary 

In this study, we performed numerical simulations with steady 
and dynamic spiral models and compared the locations of dense 
gas around spiral arms. The steady spiral models predict that the 
dust lanes are located downstream of the stellar spiral in the in¬ 
ner regions and shift upstream of the spiral in the outer regions. 
In contrast, the dynamic spiral models show no systematic ra¬ 
dial dependence of the arm-gas offsets, since the gas motions 
follow large-scale colliding flows to form spiral arms. Thus, 
our results suggest that a radial profile of arm-gas offset angle 
can be used to distinguish between the two spiral models. 

To measure the radial distributions of arm-gas offsets, 
spatially-resolved wide-field maps of old stars and cold gas are 
required. Kendall et al. (2011) used stellar mass density maps 
and gas shocks traced by the Spitzer Infrared Array Camera 
(IRAC) 8 pm data, and investigated arm-gas offset profiles of 
nearby spiral galaxies. Forthcoming integral field spectroscopic 
surveys such as the MaNGA Survey (Bundy et al. 2015) and 
wide-field mapping by ALMA will provide us accurate spatially 
resolved stellar and gas maps of nearby galaxies, and will facil¬ 
itate discrimination between the two spiral theories. 

However we should note that the steady and dynamic mod¬ 
els are not the only options. For example, Dobbs et al. (2010) 
presented a tidally-induced spiral model in which the pattern 
speed is slower than the galactic rotation everywhere but de¬ 
creases as a function of radius. In order to determine the na¬ 
ture of spiral structures in a spiral galaxy, a single test is not 
enough. Other methods, such as CO-Ha offset measurement 
(Egusa et al. 2009) and velocity-field modeling (Kuno & Nakai 
1997; Miyamoto et al. 2014), are complementary to the test pro¬ 
posed in this paper. 

The authors thank the anonymous referee for careful read¬ 
ing. IB thank also Takayuki R. Saitoh, Jin Koda, and Keiichi 
Wada for their constructive comments. Calculations, nu¬ 
merical analyses and visualization were carried out on Cray 
XT4, XC30, and computers at Center for Computational 
Astrophysics, National Astronomical Observatory of Japan. JB 
was supported by HPCI Strategic Program Field 5 “The origin 
of matter and the universe” and JSPS Grant-in-Aid for Young 
Scientists (B) Grand Number 26800099. 
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